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摘要 : 时 域 脉冲 机 制 雷达 广泛 应 用 于 地 球 探测 、 月 球 与 深 空 探测 等 领域 ， 研 究 雷 达 的 探 
测 深度 有 利于 分 析 目 标 是 否 在 雷达 探测 范围 内 。 传 统 的 计算 方法 大 多 基于 经 典 雷达 传输 方 
程 ， 在 实际 计算 时 ， 需 要 对 地 下 介质 结构 以 及 介质 特性 做 先 验 性 假设 ， 进 而 给 出 理论 的 穿 透 
深度 。 为 了 克服 传统 计算 探测 深度 方法 的 局 限 性 ， 给 出 了 一 种 新 的 计算 方法 ， 该 方法 不 需要 
对 探测 介质 做 先 验 假设 ， 直 接 从 雷达 实测 数据 出 发 ， 通 过 计算 雷达 数据 间 的 相关 性 ， 给 出 雷 
达 的 探测 深度 。 给 出 了 两 种 具体 的 计算 方法 ， 一 种 为 结合 子 波 形式 ， 另 一 种 是 利用 道 相关 形 
式 。 分 别 对 两 种 方法 进行 了 介绍 ， 并 通过 仿真 验证 了 两 种 方法 在 计算 雷达 穿 透 深度 时 的 有 效 
性 ， 同 时 分 析 了 两 种 计算 方法 存在 的 局 限 性 ， 该 研究 为 后 期 处 理 时 域 脉冲 雷达 穿 透 深度 问题 


m 提供 了 解决 方法 。 
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一 时 域 脉冲 雷达 具备 快速 、 无 损 、 高 分 辩 率 等 特点 ， 是 一 种 应 用 广泛 的 浅 层 电磁 探测 手段 。 在 地 球 
一 探测 中 ， 时 域 脉冲 雷达 广泛 应 用 于 工程 与 环境 地 球 物理 勘查 、 工 程 质量 检测 、 水 文 地 质 检测 、 基 地 考 
察 等 领域 2 ， 嫦 娥 三 号 月 球 车 搭载 的 测 月 雷达 正 是 基于 时 域 脉冲 机 制 ， 其 科学 目标 为 探测 月 壤 厚 度 
- 及 次 表层 的 结构 特性 ， 基 于 嫦娥 三 号 的 雷达 数据 处 理 也 取得 了 -一些 初步 的 结果 [+ 。 该 雷达 工作 原 
理 为 发 射 天 线 发 射 高 频 脉冲 电磁 波 ， 电 磁 波 在 传播 过 程 中 ， 由 于 地 下 介质 存在 电 性 参数 差异 ， 会 对 电 
磁 波 的 振幅 、 波 形 和 频率 等 特性 产生 影响 ， 通 过 分 析 接 收 天 线 接收 的 回 波 信号 推断 地 下 介质 结构 以 及 


介质 的 介 电 特性 。 
e 雷达 究竟 能 探测 多 深 ， 需 探测 的 地 下 目标 是 否 在 雷达 系统 的 有 效 测 距 范围 之 内 ， 在 雷达 数据 中 发 
= 现 的 特征 信息 是 否 在 雷达 有 效 数 据 范 围 之 内 ， 这 是 实际 工程 中 必须 考虑 的 一 个 重要 问题 "” 。 传 统 的 
FT 计算 方法 “" ， 大 多 需要 预 估 地 下 的 分 层 结 构 以 及 各 分 层 介质 的 物理 特性 ， 然 后 结合 经 典 的 雷达 传输 


方程 对 雷达 的 探测 深度 进行 倍 计 。 由 于 传统 探测 深度 的 计算 方法 需要 做 先 验 性 假设 ,没有 结合 雷达 的 
探测 数据 ， 所 以 只 能 给 出 一 个 理论 的 估计 值 。 本 文 提出 了 一 种 基于 相关 系数 法 雷达 极限 探测 深度 的 计 
算 方 法 ， 该 方法 不 需要 对 地 下 介质 进行 假设 ， 直 接 从 雷达 实测 数据 的 角度 出 发 ,通过 数据 间 的 相关 
性 ， 给 出 雷达 的 极限 探测 深度 。 文 中 对 于 相关 系数 法 给 出 了 两 种 具体 的 计算 方法 ,第 1 种 为 子 波 法 ， 
该 方法 利用 雷达 发 射 子 波 与 雷达 数据 进行 相关 计算 ， 得 到 相关 系数 矩阵 ， 在 相关 系数 和 矩阵 中 寻找 雷达 
的 极限 探测 次 度 ; 第 2 种 为 道 相 关 法 ， 该 方法 不 需要 子 波 ， 利 用 相 邻 道 之 间 的 相关 性 得 到 相关 系数 矩 
阵 ， 在 相关 系数 矩阵 中 寻找 雷达 的 极限 探测 深度 。 

需要 指出 ， 基 于 相关 系数 法 确定 的 极限 探测 次 度 为 雷达 实测 数据 中 最 后 一 个 有 用 回 波 信号 (强度 
大 于 噪声 ) 与 噪声 的 分 界 位 置 ， 在 该 位 置 以 下 ， 信 和 号 成 分 以 噪声 为 主 ， 即 使 存在 有 用 信号 ， 也 渡 没 在 
噪声 之 中 。 对 于 有 用 信和 号 (强度 小 于 噪声 ) 隐藏 在 噪声 中 的 情形 ， 有 用 信和 号 的 相关 性 已 经 被 噪声 破坏 ， 
采用 该 方法 无 法 提取 ， 该 情形 不 属于 本 文 的 研究 内 容 。 另 外 ， 有 些 学 者 提出 结合 有 效 的 数据 处 理 方 
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法 "中 ， 也 可 以 对 弱 信 号 、 薄 层 进行 提取 ， 但 是 很 多 时 候 需 要 信号 满足 一 些 特定 的 约束 条 件 。 
本 文 详细 介绍 了 两 种 计算 方法 的 理论 基础 及 计算 步骤 ， 并 通过 模型 仿真 验证 和 对 比 了 两 种 计算 方 
法 的 有 效 性 ， 并 对 仿真 结果 进行 了 分 析 和 讨论 。 


1 相关 系数 法 


本 节 介 绍 两 种 基于 相关 系数 对 雷达 极限 探测 深度 进行 计算 的 方法 ， 子 波 法 和 道 相关 法 。 子 波 法 在 
计算 时 需要 先 对 子 波 进行 测量 和 估计 ， 而 道 相关 法 不 需要 子 波 ， 直 接 利 用 实测 数据 。 
1.1 子 波 法 

子 波 法 计算 雷达 探测 深度 的 理论 基础 是 当 利 用 雷达 子 波 与 单 道 雷 达 数 据 进 行 滑动 相关 时 ， 真 实 的 
反射 信号 与 子 波 存在 高 相关 性 ， 骂 声 信号 与 子 波 进行 相关 运算 时 相关 性 较 低 ， 在 相关 和 矩阵 中 ， 找 到 最 
后 一 个 高 相关 信号 的 位 置 就 是 雷达 的 极限 探测 深度 。 

子 波 的 获取 可 以 用 模拟 仿真 和 实际 测量 两 种 方法 。 结 合 软件 模拟 仿真 时 ， 首 先 搭建 仿真 模型 ， 将 
发 射 天 线 与 接收 天 线 置 于 自由 空间 ， 固 定好 天 线 间 隔 ， 模 型 的 边界 采用 吸收 边界 ， 发 射 脉 冲 选 用 雷 克 


wi 子 波 ， 可 以 把 接收 天 线 接收 的 波形 作为 子 波 。 实 际 测量 时 ， 可 以 选 在 微波 暗室 或 者 空旷 的 场地 ， 尽 量 
> 保证 周边 没有 反射 介质 ， 同 样 将 接收 天 线 接收 的 波形 作为 子 波 ; 同时 也 可 以 在 获取 实测 数据 后 ， 采 用 
N 子 波 提取 法 “在 实测 雷达 数据 中 提取 子 波 。 


- 已 知 ， 子 波 wlet 长 度 为 n， 实 测 雷 达 数 据 data 共有 m 道 数据 ， 每 道 数 据 有 n. data 个 采样 点 ; 利 
一 用 子 波 法 进行 探测 深度 计算 的 具体 步骤 为 

| (1) 确 定 雷 达 数 据 中 每 个 采样 点 对 应 的 处 理 区 间 x, 
= data(1: n, j), i < (n_wlet - 1)72 

x; ; 7 4data(n. data -n + 1; n data, j), i > n data - (n - 1)72, (1) 
| data(i 一 (n -1)/2: i+ (n -1)/2, j), else 


其 中 ,i 为 第 i 个 采样 点 ; 7 为 第 j 道 数 据 ; data(1: n, j) 为 选取 第 j 道 数据 中 的 第 1 到 第 ”个 采样 点 。 
(2) 计 算 子 波 均 值 wet 以 及 每 个 采样 点 对 应 采样 区 间 的 均值 xy : 


;， 区 间 大 小 与 子 波 宽 度 一 致 ; 


J 


1 n 1 n 
wlet' = 一 >》 wlet;, x; ;' -—M x, T€ (2) 
n i=1 i 


n i=1 


ST (3) 将 子 波 与 采样 点 处 理 区 间 进 行 相关 计算 得 到 相关 系数 ，c(i, j) 为 第 j 道 数据 、 第 ;点 的 相关 系 

w 数值 : 
25 (x, ~ o, g Qulet; = wlet") 

dijs- : (3) 
IPTE 


(4) 遍 历 该 道 数据 所 有 采样 点 以 及 所 有 雷达 数据 道 ， 得 到 相关 系数 矩阵 Co 

(5) 在 相关 系数 矩阵 的 每 列 中 找到 最 后 一 个 高 相关 的 脉冲 的 位 置 ， 确 定 雷 达 极 限 探测 深度 。 
1.2 道 相 关 法 

道 相关 法 与 子 波 法 不 同 ， 在 利用 道 相关 进行 探 地 雷达 深度 估计 时 ， 不 需要 提取 雷达 子 波 ， 仪 依据 
雷达 实测 数据 ， 考 虑 到 雷达 行进 的 速度 较 慢 ， 相 邻 两 道 对 应 地 下 的 介质 结构 相似 ， 所 以 有 效 的 反射 信 
号 在 雷达 相 邻 道 之 间 具 有 较 高 的 相关 性 ， 而 噪声 为 随机 的 ， 在 相 邻 道 之 间 的 相关 性 较 低 ， 所 以 道 相 关 
法 的 理论 基础 是 利用 雷达 数据 中 有 用 反射 信号 在 雷达 相 邻 道 之 间 具 有 的 高 相关 性 以 及 噪声 在 相 邻 道 之 
间 的 弱 相 关 性 来 确定 有 用 信和 号 与 噪声 的 分 界 位 置 ， 即 可 得 到 雷达 极限 探测 次 度 。 

在 利用 道 相 关 法 进行 计算 时 ， 要 设 定 每 个 采样 点 的 处 理 区 间 ， 通 常情 况 下 ， 区 间 应 大 于 一 个 子 波 
宽度 ; 然后 与 子 波 法 相同 ， 采 用 滑动 的 计算 方法 ， 得 到 每 个 采样 点 的 相关 系数 值 。 
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道 相 关 法 的 具体 计算 步骤 如 下 : 
(1) 确定 每 个 采样 点 的 相关 区 间 ， 通 常情 况 下 ， 区 间 大 小 应 大 于 一 个 子 波 的 宽度 ; 


data(1: n, j), i < (n - 1)72 
4,7 data( n data -n + 1; n data, j), i > n data — (n - 1)72, (4) 
data(i —- (n — 1)/2: i+ (n - 1)72, j), else 


其 中 ,nn 为 区 间 大 小 ; i 为 第 i 个 采样 点 ; j 为 第 j 道 数据 ，data( 13: nm， 让 为 选取 第 7 道 数 据 中 的 第 1 到 


第 nn 个 采样 点 。 
(2) 计 算 每 个 采样 点 对 应 的 采样 区 间 的 均值 x,;: 


J 


1 n 
d eco» e (5) 
i-i 


2 (x, ; = xij) (xj, ja *i ai) 
l 


i , 
^ Im -À jj) ja — Xi js1 ) 
i=1 


ss (4) 遍历 所 有 采样 点 及 雷达 数据 道 得 到 相关 系数 矩阵 C。 
(5) 在 相关 系数 矩阵 的 每 列 中 找到 最 后 一 个 高 相关 的 脉冲 的 位 置 ， 确 定 雷达 极限 探测 深度 。 


(6) 


二 2 仿真 验证 


通过 建立 模拟 仿真 模型 ， 验 证 两 种 计算 探测 深度 方法 的 有 效 性 。 

具体 的 仿真 流程 如 下 : 首先 搭建 电磁 波 传播 的 介质 模型 ， 设 置 接收 天 线 的 最 小 可 检测 幅 值 ， 当 回 波 
言 号 小 于 该 幅 值 时 ， 以 高 斯 噪声 为 主 。 利 用 设 定 的 最 小 可 检测 幅 值 ， 可 以 预先 计算 出 雷达 的 极限 探测 深 
度 作 为 真实 值 ， 然 后 再 利用 两 种 探测 深度 的 计算 方法 从 雷达 数据 的 角度 计算 极限 探测 深度 ， 并 将 计算 结 
果 与 最 小 可 检测 幅 值 算出 的 真实 值 进 行 比 对 ， 验 证 基于 相关 系数 法 计算 结果 的 有 效 性 。 

搭建 的 仿真 模型 如 图 1， 模 型 为 多 层 介质 ， 模 型 大 小 2 m x 5 m， 其 中 最 上 方 为 空气 ,依次 向 下 包 
含 了 8 种 介质 ， 其 相对 介 电 常数 也 逐渐 增加 ,分别 为 2.9、5、10、13、17、21、25 30; 天 线 包 含 一 
个 发 射 天 线 和 两 个 接收 天 线 ( Rx_A，Rx_B) ， 接 收 天 线 位 于 发 射 天 线 的 两 侧 ， 收 发 天 线 间 距 均 为 10 em, 
天 线 距 离 第 1 层 介 质 40 em。 发 射 激励 波形 采用 雷 克 子 波 ， 中 心 频率 为 400 MHz， 时 窗 为 240 ns， 采 样 
率 为 84.8 GHz。 通 过 在 x 轴 方 向 移动 天 线 ， 获 取 在 该 模型 下 的 雷达 前 面 数据 。 

理想 的 雷达 单 道 波形 如 图 2。 理想 情况 下 ， 接 收 天 线 可 以 识别 任意 小 的 回 波 信 号 ， 在 图 2 中 ,每 
个 反射 层 位 的 反射 信号 均 被 天 线 接收 ， 各 个 反射 信号 在 单 道 波 形 上 能 清晰 分 辩 。 

为 了 更 加 接近 实际 情形 (接收 天 线 存在 最 小 可 检测 功率 /幅度 ， 当 回 波 信号 强度 小 于 该 值 时 ， 信 
号 无 法 被 检测 ) ， 设 定 4 天 线 的 最 小 可 检测 信号 的 相对 强度 为 10， 如 果 回 波 信号 相对 强度 小 于 10， 则 
该 信号 不 能 被 检测 出 来 ， 同 时 整个 时 窗 内 由 均值 为 0、 标 准 差 为 5 的 高 斯 白 噪 声 填充 ; WE B 天 线 的 
最 小 可 检测 信号 的 相对 强度 为 14， 如 果 回 波 信号 相对 强度 小 于 14 则 该 信号 不 能 被 检测 出 来 ， 同 时 整 
个 时 窗 内 由 均值 为 0、 标 准 差 为 7 的 高 斯 白 噪声 填充 。 图 3 给 出 了 改进 后 4 RRA B 天 线 接收 到 的 信 
号 ， 为 了 较 好 地 显示 弱 信 号 ， 相 对 幅 值 显 示范 围 设 定 在 -50 到 50 的 区 间 。 

从 图 3 可 以 看 出 ， 由 于 A 天 线 最 小 可 检测 相对 幅度 为 10， 由 于 第 7 和 第 8 个 回 波 信号 幅度 小 于 
10， 所 以 在 A 天 线 的 单 道 接收 波形 中 ， 只 检测 到 了 前 6 个 回 波 信号 ， 且 第 6 个 回 波 位 于 4 664 点 ; B 
天 线 最 小 可 以 检测 相对 幅度 为 14， 由 于 第 6、 第 7 和 第 8 个 回 波 信号 幅度 小 于 15， 在 B 天 线 的 单 道 接 
收 波形 中 ， 只 能 检测 到 前 5 个 回 波 信号 ， 且 第 5 个 回 波 位 于 3 493 点 。4 天 线 和 B 天 线 添加 噪声 之 后 ， 
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从 图 3 中 得 出 4 天 线 的 极限 探测 深度 为 4664 点 ; B 天线 的 极限 探测 深度 为 3 493 点 。 接 下 来 结合 
波 法 和 道 相关 法 对 雷达 数据 进行 处 理 ， 计 算 两 个 天 线 的 极限 探测 深度 ， 验 证 两 种 方法 能 否 准确 地 确定 
深度 位 置 。 
2.1 基于 子 波 法 的 计算 结 

首先 ， 子 波 的 获取 可 以 通过 搭建 自由 空间 模型 得 出 。 搭 建 的 模型 示意 图 及 获取 的 子 波 如 图 4 和 图 
5。 在 图 4 中， 发 射 天 线 与 接收 天 线 的 距离 为 10 em， 与 先前 仿真 天 线 距 离 一 致 ， 天 线 置 于 自由 空间 ， 边 
界 采用 吸收 边界 。 图 5 显示 了 获取 的 子 波 波形 ， 获 取 子 波 的 宽度 为 347 个 采样 点 。 

利用 子 波 法 得 到 相关 系数 波形 如 图 6， 取 最 后 一 个 相关 信号 的 波 谷 作为 探测 深度 位 置 ， 利用 子 波 
法 得 出 4 天 线 的 极限 探测 深度 为 4685, B 天 线 的 极限 探测 深度 为 3511。 同 时 ， 由 图 6 可 以 看 出 ， 雷 
达 在 与 子 波 相关 后 ， 各 个 层 位 信息 与 原始 回 波 信号 相 比 ， 更 加 易于 识别 ， 所 以 利用 子 波 法 进行 探测 深 
度 分 析 ， 还 能 起 到 提取 精确 反射 层 位 置 的 作用 。 
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图 1 仿真 模型 示意 图 


Fig. 1 The configuration in the simulation 
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图 2 理想 的 单 道 波形 及 局 部 放大 图 
Fig.2 The ideal signal and the partial enlarged view 
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Ideal situation 
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图 3 改进 后 的 接收 天 线 接收 波形 的 对 比 图 


Fig.3 Waveform comparison after the receiving antenna is improved 
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Fig.4 The configuration in the simulation Fig.5 The wavelet 


of acquiring wavelet 


2.2 基于 道 相 关 法 的 计算 结 

首先 估计 各 采样 点 的 处 理 区 间 大 小 ， 子 波 宽度 在 300 点 左右 ， 程 序 设 定 的 各 个 采样 点 的 处 理 区 间 为 

， 大 于 一 个 子 波 的 宽度 。 图 7 为 结合 道 相 关 法 对 探测 次 度 估计 的 计算 结果 ， 取 相关 包 络 中 间 位 置 作 
E 得 出 4 天 线 的 极限 探测 深度 为 4 677 点 ，B 天线 的 极限 探测 深度 为 3485 点 。 
2.3 结果 分 析 

相关 系数 法 ， 无 论 是 道 相 关 还 是 子 波 法 都 是 通过 相关 性 对 有 用 信和 号 进行 放大 ， 噪 声 信号 进行 压 
缩 ， 在 得 到 的 相关 系数 曲线 中 确定 最 后 一 个 有 用 信和 号 的 位 置 ， 对 比 图 3、 图 6 和 图 7， 有 用 信号 较 周 
Be A E E 这 验 
证 了 相关 系数 法 的 有 效 性 。 至 于 如 何 确 定 最 后 一 个 高 相关 信号 ， 对 于 分 层 清晰 、 回 波 信号 结构 完整 的 
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情形 ， 可 以 直接 从 相关 曲线 中 确定 ; 对 于 结构 复杂 ， 回 波 信和 号 混在 在 一 起 的 情形 ， 可 根据 具体 情况 采 
用 一 阶 偏 导 求 极 值 ， 寻 找 第 1 个 过 0 点 以 及 采用 统计 法 统计 得 出 国 值 等 方法 。 


T T T T T T 
i i ' D D 1 


Antenna A 


| 


| | | | X: 4685 i i 
-0.5 ---------]]-------]L---- — | ESEN EE ven | ——— 
| l l LE 


Correlation coefficient 
Pme] 


| 1 | I I | | 
1 000 2 000 3 000 4 000 5 000 6 000 7 000 8 000 
Sample Point 


' 
' 
' 
' 
D 
| f) 


| xosu 
Y: -0.6343 


Correlation coefficient 


1 000 2 000 3 000 4 000 5 000 6 000 7 000 8 000 
Sample Point 


图 6 基于 子 波 法 的 计算 结果 


Fig.6 The results based on the wavelet 
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图 7 基于 道 相 关 法 的 计算 结果 


Fig.7 The results based on trace correlation method 


表 1 给 出 了 4 天线 和 B RRI H ZHE B Sz Ba o PR D EAA H A 23 S YE IER th REX EET 
子 波 法 在 计算 A 天 线 的 极限 深度 时 偏差 小 于 0.45%, B 天 线 的 极限 深度 偏差 小 于 0.51%。 道 相关 法 在 
计算 A 天 线 的 极限 深度 时 偏差 小 于 0.27%，B 天 线 的 极限 深度 偏差 小 于 0. 25%。 结 果 表 明 ， 两 种 方法 
在 结合 雷达 数据 进行 探测 深度 估计 时 ， 都 能 准确 计算 出 雷达 的 极限 探测 深度 ， 误 差 小 于 1%; 道 相关 
法 略 优 于 子 波 法 ， 推 测 原因 是 因为 子 波 法 需要 子 波 输入 ， 子 波 无 法 做 到 完全 一 致 估计 ， 可 能 会 带 来 一 
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些 误差 。 另 外 ， 两 种 方法 计算 的 结果 和 真实 值 略 表 1 计算 结果 对 比 
有 偏差 ， 推断 原因 一 方面 可 能 是 因为 计算 机 仿真 Table 1 Comparison between different calculated results 
精度 ,仿真 自 带 的 不 可 消除 的 误差 ， 男 一 方面 可 4 天 线 数据 | B 天 线 数据 


能 是 因为 加 入 了 较 强 的 高 斯 白 噪声 ， 对 计算 结果 PUT "m 
产生 一 些 影响 。 

基于 相关 系数 的 计算 方法 能 较 好 地 推断 雷达 
的 极限 探测 深度 。 子 波 法 由 于 在 计算 过 程 中 需要 道 相关 法 计算 值 | 4677 | 3485 
子 波 的 输入 ， 在 实际 应 用 中 ， 子 波 不 方便 获取 或 
者 不 能 准确 获取 时 ， 对 计算 准确 性 的 影响 较 大， 同时 需要 指出 ， 子 波 法 只 适合 于 地 质 结构 较 单一 的 情 
形 ， 如 果 地 下 结构 复杂 ， 多 次 反射 三 加 会 破坏 波形 ， 此 时 利用 子 波 法 往往 不 能 得 出 一 个 较 好 的 结 
道 相 关 法 仅 依据 实测 的 雷达 数据 ， 基 于 有 用 信息 在 相 邻 道 之 间 的 高 相关 性 和 噪声 在 相 邻 道 之 间 的 随机 
相关 性 的 特点 ， 能 很 准确 地 计算 出 探测 深度 ， 由 于 该 方法 不 结合 子 波 ， 不 需要 保证 雷达 回 波 信号 的 包 
络 整体 性 ， 所 以 即使 地 下 结构 较 复杂 ， 也 能 取得 较 好 的 结果 。 


子 波 法 计算 值 4 685 3511 
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针对 传统 雷达 探测 深度 计算 方法 中 存在 先 验 性 假设 的 不 足 ， 本 文 给 出 了 两 种 基于 相关 系数 法 计算 
探测 深度 的 方法 。 当 子 波 已 知 或 比较 容易 获取 ， 同 时 地 下 分 层 结构 较 单一 ， 可 以 采用 子 波 法 进行 探测 
深度 估计 ; 当 子 波 未知 ， 地 下 结构 较 复杂 时 ， 可 以 尝试 使 用 道 相关 法 。 文 中 通过 仿真 验证 了 两 种 方法 
的 有 效 性 ， 最 后 需要 指出 本 文 只 是 提出 结合 雷达 数据 进行 探测 次 度 估计 的 一 种 思路 ， 虽 然 两 种 方法 在 
仿真 时 能 取得 较 好 的 结果 ， 但 在 处 理 雷 达 数 据 时 ， 考 虑 到 信号 波形 、 空 间 环 境 、 设 备 状态 等 ， 需 要 在 
应 用 中 加 入 更 多 的 参考 因素 进行 深入 研究 。 
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Studies of Penetrating Depth of the Time Domain Pulse Radar 


Based on Correlation Coefficient Method 
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Abstract: Radar with Time Domain Pulse Mechanism is widely used to detect the Earth, the Moon as 


well as explore deep space. The research of radar penetrating depth when it is working is useful to analyze 


whether the target is within the detection range. Traditional calculation methods are mostly based on the classic 


radar transmission equations. In actual terms, priori assumptions about the subsurface structure and dielectric 


properties of the medium have to be made, then the penetrating depth can be calculated in theory. In order to 


overcome such limitations, this paper presents a new method , which does not require any priori assumptions. It 


calculates the depth by using the correlation between data from real radar detection only. This paper introduces 


two correlation coefficient methods: one is wavelet method, and the other is trace-correlation method. Using 


simulations we have verified the validity of the two methods. We also discuss the limitations of these two 


methods in the paper. Our work should serve as an important solution for penetrating depth calculation of time 


domain pulse radar in future studies. 
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